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ABSTRACT 



Aims. We observationally investigate the relation between the photoelectric heating efficiency in photodissociation regions (PDRs) 
and the charge of polycyclic aromatic hydrocarbons (PAHs), which are considered to play a key role in photoelectric heating. 
Methods. Using PACS onboard Herschel, we observed six PDRs spanning a wide range of far-ultraviolet radiation fields (Go = 100- 
10 5 ). To measure the photoelectric heating efficiency, we obtained the intensities of the main cooling lines in these PDRs, i.e., the [O i] 
63 yum, 145 yum, and [C n] 158 /jm, as well as the far-infrared (FIR) continuum intensity. We used SpitzerflRS spectroscopic mapping 
observations to investigate the mid-infrared (MIR; 5.5-14/jm) PAH features in the same regions. We decomposed the MIR PAH 
emission into that of neutral (PAH ) and positively ionized (PAH + ) species to derive the fraction of the positively charged PAHs in 
■ each region, and compare it to the photoelectric heating efficiency. 

Results. The heating efficiency traced by ([Oi] 63 /jm + [Oi] 145 yum + [Cn] 158 fan) / TIR, where TIR is the total infrared flux, 
ranges between 0.1% and 0.9% in different sources, and the fraction of PAH + relative to (PAH°+ PAH + ) spans from (+1 1)% to 87 
(±10)%. All positions with a high PAH + fraction show a low heating efficiency, and all positions with a high heating efficiency have 
a low PAH + fraction, supporting the scenario in which a positive grain charge results in a decreased heating efficiency. Theoretical 
estimates of the photoelectric heating efficiency show a stronger dependence on the charging parameter y = GoT l/2 /n e than the 
observed efficiency reported in this study, and the discrepancy is significant at low y. The photoelectric heating efficiency on PAHs, 
traced by ([O i] 63 yum + [O i] 145 /an + [C n] 158 yum) / (PAH-band emission + [O i] 63 fim + [O i] 145 pm + [C n] 158 /an), shows a 
much better match between the observations and the theoretical estimates. 

Conclusions. The good agreement of the photoelectric heating efficiency on PAHs with a theoretical model indicates the dominant 
' contribution of PAHs to the photoelectric heating. This study demonstrates the fundamental role that PAHs have in photoelectric 

, heating. More studies of their charging behavior are crucial to understand the thermal balance of the interstellar medium. 

Key words. HII regions - ISM: lines and bands - photon-dominated region (PDR) - Infrared: ISM 



X 



1. Introduction tio of the energy emitted in the gas cooling lines ([Oi] 63//m 

. , . . ... ..... and [Ch] 158 urn are the strongest ones from Av < 5 in PDRs) 

g ; Photoelectric heating is a major heating process in photodissoci- ^ ±e ^ uy absorbed by dust grains, which is 

ation regions (PDRs), and its efficiency (e pe ) is one of the key pa- feed fe ^ dust infrared (TR) emission Previous observa _ 

rameters to understanding the energy balance there. Theoretical ^ ghow a wMe yariation of ^ sourc£s> { 

investigations suggest that small grains, in particular polycyclic from 1Q -4 in W4 9N which is illuminated by an intense uv 

aromatic hydrocarbons (PAHs), play a dominant role in pho- _ , , i., . , | h nm K » 1 oo? • <-u i *• i i n\r 

toelectric heating (Bakes & Tielens 1994) e is defined as the field < Vastel et alj | 200l|) to 1-2% m the relatively lovv-UV ex- 

toetectnc neating quaices .& iens || i^ ). e pe is aerrnea as ne PDR of the Horsehead Nebula dGoicoechea et al.l2009l). al- 

fraction of energy absorbed by dust that is converted into ki- , , ■ • , ■ • , , — : — „ T An jJ T — n 

„.f . , , . , , though mechanical heating also plays a role in W49 (Nagv et al. 

netic energy of the elected electrons and therefore into gas heat- „, * ... — — — 4 — . ^rJ?.. , , ■ r £ — n ^ 

„, aj . . , .... • 2012). Mizutam et al. 2004) showed a variation of e De from 0.06 

ing. Observationally it can be estimated by measuring the ra- — 7%Jz — 7^, — TX, j. , „ ■ » T . , i 

J J 6 to 1.2% across the 40' x 20' area of the Carina Nebula, and an 

* Herschel is an ESA space observatory with science instruments anti-correlation between e pe and the intensity of the local UV ra- 

provided by European-led Principal Investigator consortia and with im- diation field. For extragalactic sources, |Mookerjea et al.| (|201l|) 

portant participation from NASA. This work is based in part on obser- presented e pe of 0.3-1.2% toward one H n region in M 33 at a 

vations made with the Spitzer Space Telescope, which is operated by resolution of 50 pc. These variations have been attributed to dif- 

the Jet Propulsion Laboratory, California Institute of Technology under ferences in the mean charge state of the grains, i.e., a positive 
a contract with NASA. 
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grain charge results in a decreased efficiency, and the correlation 
with the intensity of the UV radiation field supports this inter- 
pretation. However, considering two properties of different dust 
populations together, large grains that are the main carriers of 
the dust IR emission and small grains that are expected to be the 
main c ontributors of the p hotoelectric effect, limits the interpre- 
tation. Habart et al. (2001) investigated the photoelectric heat- 
ing efficiency on small grains and indicated that the observed 
spatial distribution of [Oi] 63 pm and [Cn] 158 pm in L1721 
can be better explained by a mo del with a varying abundance 
of small grains across the cloud. Joblin et al] (1201 Ot) concluded 
that the evolution of PAHs and very small grains at the bor- 
der of PDRs should be considered to model the gas energetics. 
Before Herschel, estimating the gas-cooling energy in spatially 
resolved PDRs was only rarely possible because of the poor 
spatial resolution in the far-infrared (FIR) wavelength range. 
To investigate the charge state of PAHs, mid-infrared (MIR) 
spectroscopic observations are nee ded, which can be provided 
by the Infrared Spectrograph (IRS; Houck et"aT1l2004l) onboard 
the Spitzer Space Telescope dWerner et al.l 20041) . The most evi- 
dent spectral difference between neutral PAHs (PAH ) and ion- 
ized PAHs (PAH + ) is the strength of the 6-9 pm complex rel- 
ative to the 1 1.3yum feature, which is supported by l aboratory 
experi ments dDeFrees et al.lll993t iPauzat et al.lll994t iLan ghoff 
1996) and t heoretical calculation s dSzczepanski & Valall993blaF 
Hud gins & Allamandolalll995 a b). Strong variations of the ob- 
served intensity of the individual bands have been seen in 
different environments, reflecting an evolution of the charge 
state of PAHs in different physical conditions ( Galliano et al] 
2008; ISakon et al.112001 iPeeters et al.ll2002t iJoblin et al.lll996l) 7 
Recently. iJoblin et alJ (120081) and lPilleri et al] (l2012bl) proposed 
an alternative method of deriving the fraction of PAH + using a 
spect ral decomposition approach based on a few template spec- 
tra dBerne et al.ll2007t iRapacioli etai1l2005l) that include PAH , 
PAH + , and evaporating very small grains (eVSGs). 

The combination of the obser vations by Photqdete ctor Array 
Camera a nd Spectrometer (PA CS; |Poglitsch et al.l20ldl) onboard 
Herschel (Pilb ratt et al.l l2010) and the IRS onboard Spitzer en- 
ables us to investigate the relation between e pe and the fraction of 
ionized PAHs in spatially resolved PDRs. In this paper, we report 
the results for six P DRs studied in the WAD I (Warm And Dense 
Interstellar medium; lOssenkopf et al.l201 ll) guaranteed time key 
program of Herschel. WADI is aimed to investigate the physics 
and chemistry of PDRs and shocked regions with a wide range 
of physical properties. In this study, we investigate six PDRs lo- 
cated in five different regions that show clear detections of [C n] 
and [Oi] with PACS and were observed with the IRS. 



2.1.1. Horsehead 

The Horsehead Nebula emerges from the western edge of L 1630 
as a dark cloud at visible wavelengths, and its outer edge is de- 
lineated by a b right and narrow filament in the MIR emission 
(lAbergel et al.l l2003, Fig.flJ). The exciting source cr Ori illumi- 
nates the filament from the west, and a steep change in column 
density marks the western edge of the filament, while on the 
eastern side the MIR emission decreases because of the extinc- 
tion o f the incident radiation by dense material (Abergel et al. 
2003) . Using obser vations of the H2 1-0 S(l) line emission, 
lHabart et a l. (2005) proposed a model with a density gradient 
from 10 4 to 10 5 cm"" 3 with a sca le length of about 10" of this 
filament. Goicoechea et al. (2009) mapped the [Oi] 63 pm emis- 
sion at a low spectral resolution and concluded that a nonlocal 
and non-LTE treatment can be important to model the [0 1] emis- 
sion. They derived a value of e pe of 1-2%. 

2.1.2. Ced 201 

Ced 201 is a reflection nebula, illuminated by the B9 .5 star BP 
+69° 1231. Based on radial velocity measurements. IWitt et al.l 
(1 19871) suggested that this reflection nebula is the result of an 
accidental enc ounter of a small de nse molecular cloud with an 
unrelated star. Kem per et al.l (fl999) modeled cooling lines of this 
PDR and concluded that heating by PAHs and VSGs is required 
to exp lain the observed cooling-line emissions. Young lDwl et al .1 
(120021) estimated e pe to be < 0.3% based on low spatial resolu- 
tion KAO observations of [Cn] 158 pm and the upper limit of 
[0 1] 63 pm. 

2.1.3. NGC7023 

NGC 7023 is a prototype PDR, which has been widely stud- 
ied at many wavelengths. The UV radiation from the exciting 
B star, HD200775, creates three main PDRs in this nebula; the 
brightest NGC 7023 North- West (hereafter NGC 7023 NW) 
is located about 40" northwest of the star, another PDR lies 
about 70" south (NGC 7023 S), and NGC 7023 Ea st (hereafter 
NGC 7 023 E) is loc a ted abo ut 170" east of the star dBerne et al.l 
120071) . iPilleri et all (l2012bl) showed a clear difference in the 
spatial distribution of PAH and PAH + in these three PDRs. 
IJoblin et al] d2010l) presented the first results from HIFI observa- 
tions along a cut through NGC 7023 NW and S, suggesting that 
both the [Cn] emission and the aromatic infrared-band (AIB) 
emissions in the MIR arise from the regions located in the tran- 
sition zone between atomic and molecular gas, providing new 
insights into the importance of the PAH charge evolution in the 
energetic studies of PDRs. 



2. Observations and data reduction 

2.1. Targets 

Table Q] shows the global properties of the investigated re- 
gions. They cover three orders of magnitude in far-ultraviolet 
(FUV; hv = 6-13.6 eV) flux, Gq in units of the Habing field 
(1.6 x 10" 6 Wnr 2 lHabinglll968l) . Figure Q] shows the PACS 
field-of-view and the extracted area for each PDR. Some of 
these sources have been studied and modeled in detail based 
on the observations with the Herschel-Heterody ne Instrument 
for the Far-Infrare d (HIFL lde Graauw et al.ll201of) in WADI (see 
Joblin et al l I20T0I for NGC 7023 and iFuente etall I20T0I and 
Pille rietal.ll2012al for Mon R2). 



2.1.4. Carina 

The Carina Nebula is a massive star-forming re gion complex 
with 65 O-type stars at a distance of 2.3 kpc (Smith 2006), 
which provides the nearest example of a very massive star- 
forming region and ha s been observed at many wavelengths 
(ISmith & Brooksl2007l) . Trumpler 14 and 16 are the most promi- 
nent clusters. The PDR properties have b een investigated us- 
ing FIR and submil i meter emission lin es dKramer et al.l [2008; 
Mizutani et al. 1 12004 iBrooks et alj|2003l) . The region excited by 
Trumpler 14 ha s a FUV flux of G o = 7 x 1 3 and a high de nsity of 
2 x 10 5 cnT 3 dKramer et al.ll2008l) . lPreibisch et al.l (1201 ll) present 
a large deep 870 pm continuum map, showin g that the total mass 
is about 2 x 10 5 Af™ . iGaczkowski et all (120 1 31) recently published 
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Fig. 1. PACS footprints of the blue camera ([O i] 63 fim) and the red camera ([O i] 145 fim and [C n] 158 fim) overlaid on IRAC 8 fim, 
except for (f) Mon R2. (There IRAC 8 fim is saturated, and the total intensity of AIB bands derived from fitting of IRS spectra (see 
text) is shown.) Green circles or boxes show the areas from which our PAC S and IRS spectra a re extracted (subregions in Table [3) . 
For Mon R2 (f), the positions of four IR sources (IRS 1—4) and PDR 1 from lBerne et al.1 ( 2009) are marked. 
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Table 1. General properties of the observed targets. Go only gives typical values within the individual regions. See Table|4]for local 
variations. 



Object 




Distance 


Exciting 


source 




Reference 






[kpc] 


Name 


Spectral Type 


[Habing field] 




Horsehead 




0.41 


crOri 


09.5V 


100 


1,2 


Ced 201 




0.42 


BD+69° 1231 


B9.5V 


300 


3 


NGC 7023 


E 


0.43 


HD 200775 


B3Ve-B5 


250 


4 




NW 


0.43 


HD 200775 


B3Ve-B5 


2600 


4 


Carina 


N 


2.35 


Trumpler 14 


03 


7000 


5 


Mon R2 




0.85 


IRS1 


B0 


5 x 10 5 


6,7 



a See text (Sect.l2.1[> for the defi n ition. 

Reference: ' iMenten et all d2007h . 2 lHabart et alJd2005h . 3 lYoung Owl et alJ < l2002h . 4 lPilleri et al.lll2012bh . 5 lKramer et aIU2008h . 6 ll)ownes et al] 
l fl975T ), 7 iRizzo et al] imm 



PACS and SPIRE continuum maps of the Carina Nebula com- 
plex to study embedded young stellar objects. In this study, we 
analyzed a PACS position in Carina North (hereafter Carina N), 
excited by Trumpler 14. 

2.1.5. Mon R2 

Mon R2 is a close-by ultracompact H n (UCH n) region at a dis- 
tance of 850 pc, which comprises several PDRs that can be spa- 
tially resolved at both mm and IR wavelengths. The most intense 
UV source, which coincides with the infrared source IRS1, is lo- 
cated at the center of the cometary shaped UCH n. The brightest 
PDR of Mon R 2 is illuminated by an extremely intense UV field 
(Go ~ 5 x 10 s : IRizzo et alj f2003). and its PAH emission peaks 
at about 20" northwest of IRS1 (Fig.Q}. Several other PDRs lie 
around IRS 1 and span different physical conditions such as tem- 
perature, density, column density, and UV f ield (e.g. lRizzo et al.l 
l200llBerne et al.l2009HPilleri et al.l2012ab . All PDRs have very 
bright MIR spectra, consisting of the emission from AIBs, H2 
rotational lines, and the continuum. The spatial extent of these 
PDRs (~ 0.03 pc) yields an angul ar size of ~ 8", a s shown by 
PAH and H 2 emissions in the MIR dBerne et al.l2009h . The shape 
of the MIR spectra varies at different PDRs, reflecting the photo- 
processing of the AIB carriers. 

2.2. Far-infrared spectroscopy with Herschel 

The FIR spectroscopic observations were performed with PACS 
as part of the WADI program. The observational parameters are 
listed in Table [2] All measurements are single pointed observa- 
tions with the PACS spectrometer consisting of a 5 x 5 array of 
spatial pixels (spaxels) with a size of 9.4" each (see Fig. |T}. We 
pipelined the data from leve l to 2 using the Herschel data pro- 
cessing system (HIPE: [Otil |20 1 Oh version 8. We converted the 
unit of the pipelined spectra, Jy/spaxel, to Wm^yUirT 1 sr~' us- 
ing the field of view of a spaxel (9.4" x 9.4"). For the Chop/Nod 
observations, we checked the contamination at the OFF positions 
for [0 1] 63 /im, 145 yum, and [C n] 158 /im using the script pro- 
vided in HIPE (ChopNodSplitOnOff.py), which provides sep- 
arate spectra for the ON and OFF positions. Only the [Cn] 
158^m emission in the Horsehead shows significant OFF con- 
tamination. Since the northern side of the OFF positions shows 
stronger [Cn] 158/mi, and the emission from the southern side 
can be attributed to a diffuse component that is not associated 
with the region, we used only the OFF measurement of the 
southern side to obtain the flux from the source. For the un- 
chopped grating-scan observations, some emission is detected 



at the OFF positions in [Cn] 158yum (Carina N and Mon R2) 
and [Oi] 63 fim (Carina N). The line intensities at the OFF po- 
sitions are 4-7%, 2-5%, and \-\% of that at the ON positions. 
We subtracted the OFF emission in Mon R2 but not in Carina N, 
because the detection of [0 1] 63 /im in the Carina N OFF po- 
sition indicates a real contamination from local dense clouds, 
while the [Cn] 158 fim emission traces more diffuse region and 
can be attributed to the large-scale diffuse Galactic [C n] emis- 
sion. To determine the continuum levels, we subtracted the OFF 
measurements in all observations to remove the telescope back- 
ground. 

The uncertainty of the absolute flux calibration is 11-12% 
for the line and continuum emission^. A larger uncertainty is 
expected in weak sources for unchopped grating scans. The 
reproductivity in total absolute flux (telescope + source) of 
the unchopped grating scan is 4% (peak-to-peak). The rela- 
tive uncertainty against the source flux can be expressed as 
0.04 x ^(telescope + source)/F(source), where F is the flux. For 
Mon R2, this means an uncertainty of < 13% and < 5% for the 
continuum at 105-180jt/m and at the peak of [Cn] 158 fim. For 
Carina N, it is < 10% and < 5% at the peak of the [Oi] 63//m 
and [C 11] 158 jum emission lines and 20-50% for the underlying 
continuum level. 

We defined the region to be used in our study depending on 
the morphology of each source so that the widest range of phys- 
ical conditions is included. Since our PACS observations are not 
fully sampled, we cannot exploit the entire spatial information. 
Instead, we selected a few typical areas in each source and ex- 
tracted the [0 1] and [C 11] line intensities and the continuum flux 
from PACS observations and combined them with the IRS re- 
sults. These areas are listed in Table [3] shown as green boxes or 
circles in Fig.Q] and are explained in detail in the following. We 
extracted spectra based on the geometrical area and did not apply 
a beam-size correction. The FWHM of the PACS spectrometer 
beam at 63//m and 158//m is ~ 9" and ~ 11.5", respectively, 
and the uncertainty from the difference of the beam size is esti- 
mated to be ~ 15% (see Appendix lAb. 

For Horsehead, Carina N, and Mon R2, all observations were 
made with almost the same position angle, i.e., the PACS spax- 
els observed almost the same area in the sky for different wave- 
lengths (Fig. [TJ. In this case, we extracted the areas to be ex- 
amined on the basis of the PACS spaxels. In the Horsehead, the 
spaxels are aligned to the ridge. We defined three areas; the first 
one lies at the western side of the ridge in the ionized gas (an av- 
erage over four spaxels except for the north-end spaxel), the sec- 
ond area is along the ridge (over five spaxels), and the third area 

1 PACS spectroscopy performance and calibration document v2.4 
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Table 2. PACS observation summary. 



Ohiert 




RA fROOO"! 




Obs. ID 


fine morlp 




Horsehead 




5 h 40 m 53 s .7 


-2°28'00".0 


1342228244 


LineSpec, Chop/Nod 


[Oi] 63 /mi 










1342228245 


LineSpec, Chop/Nod 


[Cn] 158 yum 












T " O ^"11 /XT J 

LineSpec, Chop/Nod 


[Oi] 145 /mi 


Ced 201 




22 h 13 m 24 s .6 


+70°15'04".6 


1342211847 


LineSpec, Chop/Nod 


[Oi] 63 /mi 










1 i /no i tnno 


T * o /"il . /XT 1 

LineSpec, Chop/Nod 


LOi] 145 /mi, LLnJ 158/mi 


NGC 7023 


E 


21 02 09 .4 


+68 09 15 .3 


1342209729 


T " O /"II /XT J 

LineSpec, Chop/Nod 


[Oi] 63 /mi 










134zZZzZ31 


LineSpec, Chop/Nod 


LuiJ 143 /jm, LCiiJ Ijs/jm 




NW 


21 h 01 m 32 s .4 


+68°10'25".0 


1342208910 


LineSpec, Chop/Nod 


[Oi] 63 /mi 










1342222230 


RangeSpec, Chop/Nod 


[Oi] 145 /mi, [Cn] 158/im 




NW 


21 h 01 m 35 s .7 


+68°10'01".0 


1342196683 


LineSpec, Chop/Nod 


[Oi] 63 /mi 










1342197033 


RangeSpec, Chop/Nod 


[Oi] 145 /mi, [Cn] 158/mi 


Carina 


N 


10 h 43 m 25 s .7 


-59°34'35".3 


1342214631 


LineSpec, Unchopped 


[Oi] 63 /iin, [Cn] 158 /mi 


Mon R2 




6 h 07 m 46 s .2 


-6°23'08".3 


1342228453 


LineSpec, Chop/Nod 


[Oi] 63 /mi 










1342228454/6 


RangeSpec, Unchopped 


[Oi] 145 yt/m, [Cn] 158/mi 



Table 3. Line intensities observed by PACS. Subregions are shown by green circles or boxes in Fig. [TJ and are described in detail in 
Sect.O 



Object Subregion 


Line intensities [Wm 1 sr '] 
[Oi]63yum [Oi] 145pm [Cn] 158/mi 


Horsehead 1 H n region side 

2 along the ridge 

3 molecular side 


(1.01 ±0.02) x 10-' (1.28 ± 0.03) x 10" 8 (9.49 ± 0.23) x 10^ 
(2.39 ± 0.02) x 10- 7 (5.08 ± 0.03) x 10~ 8 (2.27 ± 0.02) x 10~ 7 
(1.29 ± 0.02) x 10^ 7 (2.75 ± 0.03) x 10^ 8 (1.66 ± 0.02) x 10^ 7 


Ced 201 1 9.4" diameter 
2 28.2" diameter 


(6.02 ± 0.24) x 10~ s (5.48 ± 0.47) x 10"* (7.55 ± 0.09) x 10~ 8 
(4.59 ± 0.13) x 10^ 8 (4.33 ± 0.24) x 10" 9 (6.09 ± 0.04) x 10^ 8 


NGC7023 E 1 9.4" diameter 
2 28.2" diameter 


(1.53 ± 0.03) x 10-' (4.48 ± 0.07) x 10~ 8 (3.23 ± 0.03) x 10"' 
(1.56 ± 0.02) x 10^ 7 (3.90 ± 0.03) x 10^ 8 (3.04 ± 0.02) x 10^ 7 


NGC7023 NW 1 cavity 

2 interface 

3 molecular region 


(6.90 ± 0.06) x 10-' (6.66 ± 0.21) x 10- 8 (5.53 ± 0.11) X 10"' 
(2.3 1 ± 0.02) x 10~ 6 (3.73 ± 0.09) x 10~ 7 (9.09 ± 0.24) x 10~ 7 
(1.47 ± 0.01) x 10- 6 (3.05 ± 0.05) x 10- 7 (6.59 ± 0.10) x lO" 7 


Carina N 1 8 /mi clump 

2 east-west ridge 

3 southern 5x3 spaxels 


(4.89 ± 0.04) x 10-" - (4.46 ± 0.04) x lO" 6 
(4.00 ± 0.02) x 10~ 6 - (4. 10 ± 0.03) x 10~ 6 
(3.51 ±0.02) x 10- 6 - (3.98 ± 0.03) x 10^ 6 


Mon R2 1 around PDR1 

2 southwest 

3 northeast 


(7.68 ± 0.08) x 10^ 6 (1.66 ± 0.06) x 10 -6 (1.64 ± 0.07) x 10 -6 
(9.37 ± 0.12) x 10^ 6 (3.56 ± 0.14) x 10^ 6 (3.06 ± 0.12) x 10^ 6 
(9.05 ± 0.17) x 10^ 6 (2.38 ± 0.09) x lO" 6 (2.73 ± 0.10) x 10^ 6 



is in molecular gas, at the eastern side of the ridge (five spaxels). 
In Carina N, we selected three regions; a single spaxel toward a 
clump seen in the IRAC 8 /mi map, a ridge-like east-west struc- 
ture in the 8 /mi map including this clump (five spaxels), and the 
southern 5x3 spaxels where the PAH emissions in the IRS spec- 
tra are prominent. In Mon R2, we selected three spa xels around 
PDR1 (following the definition in lBerne et al. 2009), where the 
PAH + fraction is suggested to be low, and toward the southwest 
and northeast inne r edges of the PAH + distributions shown in 
Bern e et al] (|2009). They correspond to the three green boxes in 
Fig. QJ from southwest to northeast. 



For Ced 201, NGC 7023 E, and NGC 7023 NW, the [Oi] 
63 /mi line was observed in different seasons of the year than the 
[Oi] 145 /mi and [Cn] 158 /mi observations, and the difference 
of the position angle is significant. Therefore, we defined cir- 
cles as areas to be studied, and took a weighted mean of PACS 
spectra based on the geometrical overlap between the defined 
circles and the PACS spaxels. For Ced 201 and NGC 7023 E, 
two circles were defined with a common center for both PACS 
footprints and diameters of 9.4" and 28.2". For NGC 7023 NW, 
we selected three typical regions that trace diff erent ionizatio n 
fractions of PAHs based on the analysis in iPilleri et al.l d2012bl) . 
The areas were defined by circles with 12" diameter, at the in- 



terface with strong PAH emission, toward the molecular cloud, 
and in the cavity toward the exciting star (Fig. [TJ. 

After extracting the spectra in each region as described 
above, we obtained the [Oi] 63 /mi and [Cn] 158 /mi intensi- 
ties, as well as [Oi] 145 /mi when available, by a Gaussian fit 
after linear baseline subtraction (Fig. [2] Table [3}- In general, 
the rat io of [O i] 63 um / [Cn] 158/mi traces the density of the 
PDRs jRollig et alj|2006l) . Among our PDRs, NGC 7023 E and 
Ced 201 show a low ratio (< 1), indicating a low density, and 
Mon R2 has the highest ratio (> 4). This trend is roughly con- 
sistent with the density estimate in previous studies (Table [4] 
Appendix IB- The [Oi] 145/mi/63/mi is > 0.09 for all tar- 
gets and 0.3 at maximum, which exceeds the PD R model pre- 
diction for the relevant phy sical parameter ranges (Rolli get al.l 
2006; Kaufm an et al1 fT999). Possible reasons are high optical 
depths of the [O i] 63 /mi emission, a suprathermal population of 
[Oi] b y collision with H2 , and foreground absorption in [Oi] 
63 /mi (Liseau et al.l l2006h . Although the optical depth of the 
[0 1] 63 /zm line is taken into account by the PDR models for 
simple geometries, the overlapped several PDR clumps along the 
line-of-sight or the edge-on geo metry can cause more significant 
self-a bsorption of [0 1] 63 /mi (Habart et al. 2003; Okad a etafl 
120031) . and the non -local calculation suggests a h igher [Oi] 
145 /mi/63 /mi ratio (Eli tzur & Asensio Ramos 2006). The fore- 
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Fig. 2. Examples of [O i] 63 /mi, 145 /mi, and [C n] 158 /mi spec- 
tra in each PDR. The numbers labeled in figures after the name 
of the object show subregions (see TableO. 



ground absorption cannot be quantified in our PDRs because of 
the lack of velocity-resolved [O i] 63 /zm observations. However, 
the HIFI observations of the velocity-resolved [C n] 158 /mi im- 
ply only minor effects toward our PDRs. 

2.3. Mid-infrared spectroscopy with Spitzer 

We analyzed the MIR spectra observed with the Short-Low 
(SL1 and SL2) module of IRS onboard Spitzer except for 
NGC 7023 E, for which we used the I SOCAM highly-processed 
data product (Boulanger et all 120051) . We used the same data 
cubes as in Pill eri et aLld2012bl) for the Horsehead, Ced 201, and 



NGC 7023, and the data reduction is described in that paper as 
w ell. For Mon R2, the IRS data and their reduction are shown 
in lBerne et al.l (2009). For Carina N, we collected the IRS obser- 
vations fr om the Spitzer dat a archive and analyzed them using 
CUBISM dSmith et alj|2007h . To estimate the continuum, Short- 
High (SH) and Long-Low (LL) spectra were also reduced and 
extracted when available. Their intensities are scaled to match 
the SL spectra, then the whole spectra are scaled to match the 
photometric data; IRAC 8 /mi and MIPS 24/zm for Horsehead, 
Ced 201, and NGC7023 E, and IRAC 8 /mi for NGC7023 NW 
and Carina N. For Mon R2, the absolute flux scaling was not 
applied because IRAC 8 //m is saturated, and the available spec- 
tral range does not cover the MIPS 24 //m band. The correction 
should be less significant in Mon R2 because it is the brightest 
source in our samples. Applying this correction does not affect 
any trend in the following results and changes none of our con- 
clusions. We extracted a spectrum from the same area as for the 
PACS spectra described in the previous subsection. 



3. Analysis 

3.1. Total infrared intensity 

Observationally, e pe is measured as ([Oi] + JCul) / TIR, 
where TIR is the total infrared flux (3-1 100 /mi; iDale & Heloul 
2002), representing the fraction of the input energy that is con- 
verted into the coo ling lines. The total FIR flux, 42. 5-122. 5/mi 
dHelouetalJll988h or the integrated flux of the thermal emis- 
sion by large dust grains, is often used as a tracer of the input 
energy. Here we used the TIR instead, which is the emission 
from all dust grains, including PAHs and VSGs, to represent 
the total input energy. As dominant cooling lines, we summed 
the intensities of [Oi] 63 /mi, 145 /mi, and [Cn] 158 /zm, except 
for Carina N, where the [O i] 145 /zm was not observed, and we 
used ([Oi] 63 /mi + [Cn] 158 /zm). A second PACS observation 
1.3' away in Carina N, which is not used in this study because 
it does not overlap with the IRS observations, shows an [Oi] 
145/mi/63/mi intensity ratio of 0.05-0.15. If our position has 
the same ratio, neglecting the [Oi] 145 /mi emission underesti- 
mates e pe by 8% at most. 

To estimate TIR as the energy input tracer we need to inte- 
grate the grain emission of 3-1 100 /zm. The contribution from 
3-5.5 /zm is negligible because of the low intensity at 5.5 /zm in 
all regions. To obtain the energy in the range of 5.5 -14 /zm, we 
fit the I RS spectra using the procedure described in Piller i et al.l 
(2012b) and integrated the resulting fit (see Sect. 13.21 1. 

For A > 14/zm, the available data and the quality are not 
uniform, and we fine-tuned approaches from region to region 
to estimate the integrated flux (see Fig. [3), although the basic 
approach was the same; we fit the FIR spectra from PACS ob- 
servations with the thermal dust model, extrapolated the fit in 
the MIR to longer wavelengths, and connected the two fits at the 
intersection point. 

In Horsehead, Ced 201, and NGC 7023 E, MIR spectra are 
available up to 35 /zm. We estimated the integrated flux for 14— 
28 /zm by direct integration excluding the strong emission lines 
for Ced 201 and NGC 7023 E, and by a linear fit to the spectra in 
the Horsehead because of the lower signal-to-noise ratio (S/N). 
The 28-32 /zm range is fitted by a linear function, and the results 
are extrapolated to longer wavelengths. In NGC 7023 NW and 
Carina N, IRS/SH data are available for 14-19/zm, which are 
directly integrated, and we assumed a flat spectrum from 19 /zm 
up to the wavelength where the big-grain emission, peaking in 
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Fig. 3. Continuum fit at the same positions as in Fig. [2] The black lines are the observed spectra, while the red and blue lines are 
the models used to estimate the TIR. The blue asterisks are the data points that are used for fitting the FIR thermal dust emission 
(see text). The blue curve in the FIR is the fit with the thermal dust model, that in the MIR indicates either the direct integration, 
the linear fit, or the assumed constant value. The red curve below 14/mi shows the fit described in Sect. 13.21 The black asterisk in 
NGC7023 NW is the MIPS 24 yum flux. 



the FIR, exceeds this level. For Mon R2, we assumed that the 
flux at A > 14yum is constant (Fig. [3j. 

For the PACS Chop/Nod LineSpec observations (Table 0, 
the continuum levels are available around the emission lines. 
They are indicated as blue asterisks in Fig. [3] In NGC 7023 NW 
and Mon R2, the full spectral data from the PACS RangeSpec 
observations at > 70 //m, after excluding strong emission lines, 
were used to fit the FIR continuum. In Carina N, the un- 
certainty of the absolute continuum flux is too large (see 
Sect. 12.2b . Therefore we obtained the photometry data at 
blue (70/im) and red (160/itn) bands from the Herschel data 
archive. We fit the FIR conti nuum using the dust emissivity of 
lOssenkopf & Henningl (IT994). for a gas density of 10 6 cm" 3 and 
grains with thin ice mantles, considering the temperature and the 
column density as free parameters. Since we discuss only the in- 
tegrated IR flux, the choice of the dust emissivity is not critical. 
Where this fit was exceeded by the extrapolation from the MIR 
spectrum, we connected these two at the intersection point. The 
final fits are shown as blue lines in Fig. [3] 

Since the PACS observations are not fully sampled and the 
IRS spectra were extracted just using the geometrical area, the 
difference of the PSF could affect the resulting fit of the spectral 
energy distributions (SEDs). However, the uncertainty decreases 
for larger regions (average over several spaxels) and it does not 
change the general trend of e pe (Sect.|H Fig. [5j. 



In NGC 7023 NW, the MIPS 24jt/m flux is available but was 
not used to scale the IRS spectra because the spectral range of the 
IRS does not cover the MIPS 24 fim band. Figure [3] shows that 
the MIPS 24 /im flux matches the assumed constant continuum 
flux. 

We translated the obtained TIR into the impinging FUV flux. 
By assuming that the entire FUV energy is absorbed by dust 
grains and reradiated as infrared flux, we estimate Go (TIR) as 
4n x TIR/(1.6 x 10~ 6 ) (Table ID . The systematic uncertainty 
of TIR is ~ 15%, which comes from difference of the beam 
sizes at different wavelengths (Sect. l2~2l . In this formula we as- 
sumed a spherical geometry. When we resolve an edge-on PDR 
with a shorter thickness compare d to the length of th e line-of- 
sight, Go(TIR) is overestimated dMeixner et al.llT992T) . Another 
uncertainty of Go(TIR) is the contribution of photons with an 
energy outside of 6 eV < hv < 13.6 eV to the thermal heating 
of the dust. For later B-type stars, the contribution of < 6 eV 
becomes significant, wh ich causes the overestima te of Go(TIR). 
The SED of a B9 star bv lCastelli & Kuruczl d2004 indicates that 
only 1/4-1/5 of the total energy is emitted in the range of 6 eV 
< hv < 13.6 eV from the star. We did not apply the correc- 
tion and treated Go(TIR) as an upper limit. Detailed compar- 
isons of Go(TIR) with previous studies and other diagnostics are 
described in Appendix iBl and the final Go estimates, which are 
used in the following analysis, are also listed in Table [4] 
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Table 4. TIR and Go(TIR) and the parameters used to estimate the charging parameter y = GoT l ^ 2 /n e . The uncertainty of TIR and 
Gn(TIR) is ~ 15% (see text). 
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Fig. 4. Fit of the MIR spectra using the PAHTAT procedure dPilleri et alj|201 2b) at the same positions as in Fig. [2] assuming the 
mixed extinction model and Ry = 3.1. The black line shows the observed spectrum, the red line is the fitted spectrum, the blue line 
represents the continuum, green, orange, light-green and light-blue show the PAH , PAH + , PAH*, and eVSG components. 



3.2. Ionization of PAHs 

To derive the contribution of PAH + in each region, we applied 
the PAHTAT (P AH Toulouse Astron omical Templates) proce- 
dure described in IPilleri et al.l(l2012bl) . This procedure fits MIR 
spectra using a minimal set of template spectra. The PAH-related 
templates comprise PAH , PAH + , larger ionized PAHs (named 



PAH*, see lJoblin et al.ll20 08). and evaporating very small grains 
(eVSGs). In the MIR the PAH and PAH + templates are char- 
acterized by a strong difference in the relative strength of their 
7.7 yum and 1 1.3 jum band features; PAH + has a stronger 7.7 fim 
band. The PAH* population consists of large ionized PAHs, 
which was introduced to provide a better fit to planetary neb- 
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Table 5. The fraction of each component relative to the total band emission in the PAH fit. The variation in each column describes 
the range covered by the different assumptions in the fit. 
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ula spectra: its MIR spectrum is similar t o PAH + , but the 
7.7 fjm band is shif ted to longer wavelengths dJoblin et al.f 2008; 
iPilleri et alJl2012bT) . eVSGs are an intermediate population be- 
tween PAHs and classical VSGs. They present both a broad band 
and continuum emission in the MIR range. The PAHTAT proce- 
dure also allows simultaneous fitting of the gaseous emission 
lines, underlying continuum, and the extinction by dust grains 
along th e line-of-sight. More details of the fitting tool can be 
found in Piller i et all d2012b). We assumed a linear continuum 
except for Mon R2, where a continuum of two slopes connected 
at 10 pun was adopted since big grains contribute significantly to 
the MIR continuum because of the strong UV field. To explore 
the parameter space of the fitting, we used all combinations of 
the following different assumptions: the dust extinction curves 
with Ry = 3.1 and 5.5, geometries in which the absorbing ma- 
terials are placed in the foreground or are fully mixed with the 
emitting materials, and since the presence of PAH A in PDRs with 
an intermediate UV radiation field strength is doubtful, we also 
tested the results including or excluding the PAH 1 component 
from the fit. The difference between these different assumptions 
is included into the uncertainty of the results, although they are 
typically smaller than the errors estimated by the discrepancy be- 
tween the observed spectra and the model (see below). The fitted 
spectra are shown in Fig.|4] and the derived fraction of each PAH 
and eVSG contribution to the total band flux is shown in Table|5] 

In the following we define the fraction of positively 
ionized PAHs by the integrated intensity ratio as /(+) = 
f ^pah + /( § ^pah° + J/pah+X where the integration is made over 
the wavelength from 5.5 /mi to 14 /mi. The uncertainty of this 
fraction can be attributed to two different effects: ( 1 ) fitting errors 
and (2) systematic errors in the assumptions within PAHTAT, 
i.e., due to intrinsic uncertainties in the templates. To estimate 
the first contribution, we integrated the absolute discrepancy of 
the observed spectra from the model. We assumed that this resid- 
ual consists of the uncertainty of the band emissions of small 
grains (PAH , PAH + , PAH*, and eVSGs) and the continuum 
emissions; J (residual) = A(band) + A(continuum), and the ra- 
tio of the uncertainty A(band)/A(continuum) is the same as the 
ratio of the integrated flux J /band/ f /continuum- Thus, the uncer- 
tainty of the integrated band emission is defined as A(band) 
j (residual) x J 4and/(J 4and + / /continuum)- Then we attribute it 
to the uncertainty of the PAH + flux, i.e., the final uncertainty of 



/(+) is expressed as A(band)/( J I PAH ° + J IpAH + )- This is conser- 
vative, because it assumes that the entire uncertainty in the band 
emission comes from that of PAH + . 

The second uncertainty in the obtained /(+) comes from the 
templates for PAH and PAH + . The templates are constructed 
based on observational data and a mathe matical blind signal sep- 
aration (BSS) method (IPilleri et al.ll2012.bb . The spectral proper- 
ties of the BSS-extracted spectra, templates, a nd their assign- 
ment t o PAH + and PAH° have been discussed in Rapacioli et al 



20051). iBerne et all (120071). iJoblin et al J (120081) . and lBerne et al 



20091) . Recently. iRosenberg et all d201 ll) compared the BSS ex 



traction to the theoretical spectra with the density functional the- 
ory and showed good agreement, although this is for a spectral 
range of 10-19.5 fim. A precise determination of the uncertainty 
is difficult, however, given the nature of the extraction method 
and our limited knowledge of the PAH populations that e xist in 
space. Using a Monte Carlo approach, iRosenberg et al.l (1201 ll) 
showed in their BSS extraction a l<x uncertainty of typically 
10% and up to 30% in some parts of the spectrum (see their 
Fig. 3). However, this uncertainty in the templates systematically 
propagates to the uncertainty of /(+) and does not change the 
trend and our conclusion discussed below. Therefore, we present 
our errorbars in the following figures based on the fitting error 
described above. 



4. Results and discussion 

4.1. Relation between e pe and /(+) 

Figure [5] shows the relation between the e pe = ([Oi] 63/mi + 
[O i] 145/mi + [C n] 158 /an) / TIR and the fraction of positively 
ionized PAHs /(+). Both e pe and /(+) strongly vary among the 
regions; e pe varies from 0.1% to 0,9%, and /(+) from (+11)% 
to 87 (+10)%. All positions with a high /(+) show a low e pe , 
and all positions with a high e pe show a low /(+). This trend 
supports the theoretical expectation in which e pe decreases when 
grains are positively ionized, because the energy required to eject 
electrons from positively charged grains is higher than that from 
neutral grains. Here we take the ionization of PAHs as a gross 
indicator of the positive charging of grains in general, which is 
generally quantified by the charging parameter (y) as shown in 
Sect. 14.21 Figure [5] directly compares of the observationally de- 
rived e pe and /(+), independent of the use of PDR models to 
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Fig. 5. e pe = ([Oi] 63jum + [Oi] 145 jum + [Cn] 158/^m) / TIR Fig. 6. e pe = ([Oi] 63 + [Oi] 145 + [Cn] 158 pm) / TIR 



versus the fraction of ionized PAHs (/(+)). 



quantify the physical conditions, over a wide range of the phys- 
ical properties of PDRs. 

Using ([Oi] 63 pm + [Oi] 145 pm + [Cn] 158/itn) / TIR as 
a tracer of e pe contains several assumptions. We neglect the heat- 
ing by the collisional de-excitation of vibrationally excited H2. 
The H2 de-excitation heating becomes a do minant heating pro - 
cess in a dense PDR with a low UV field dRollig et al] l2006). 
Among our PDRs, only Go and the upper value of the density 
of the Horsehead are close to the regime where the contribution 
from the photoelectric heating and H2 heating is comparable. 
Those in other PDRs indicate that the photoelectric heating is 
dominant. We also neglect the cooling by H2 emission. This is 
partly justified because the H2 emissions as a consequence of 
the excitation by the UV-pumping and the H2 formation do not 
contribute to the est imate of the photoelectric heating efficiency. 
Haba rt etaD (1201 ll) showed intense H2 lines in the Horsehead 
(0-0 S(0) to S(3) and 1-0 S(l)), which are, as a sum, com- 
parable with the ([Oi] + [Cn]) intensity. In NGC 7023 E and 
Ced 201, the sum of th e H2 pure rotational emissions S(0) to 
S(3) (lHabartetal.l201ll for NGC 7023 E, and from the IRS spec- 
tra in the Spitzer archive for Ced 201) is ~ 80% and ~ 60% of 
the ([Oi] + [Cn]) intensity, respectively. Neglecting these H2 
emissions may underestimate e pe by a factor of < 2. In other 
regions, the H2 pure rotational emissions from the IRS spectra 
in the Spitzer archive show the intensity of < 20% of ([Oi] + 
[Cn]). 

As mentioned in Sect. 12.21 the [0 1] 63 pm emission is indi- 
cated to be optically thick in most targets. Therefore, the optical- 
depth-corrected [0 1] 63 pm intensity is higher than the observed 
values. However, we observe the emission that leaves the cloud, 
and this contributes to the net cooling of the cloud. This justifies 
the direct use of the observed ([Oi] + [Cn]) to express the gas 
cooling, except if the foreground absorption is significant, which 
is unlikely for our PDRs, as mentioned above. 

As discussed for Go(TIR) (Sect. 13. il l, photons with an energy 
of < 6 eV contribute to TIR, but not to the photoelectric heat- 
ing. Therefore, estimating e pe with the observed TIR provides a 
lower limit. For the exciting source(s) of spectral types earlier 
than early B, the effect is not significant in terms of the abso- 
lute factor and the variation between different spectral types. For 
Ced 201, which has an exciting source of B9.5, we may under- 
estimate e pe by a factor of a few. This uncertainty, however, does 
not change the overall trend in Fig. [5] and our conclusions. 



versus the charging parameter (y = GoT 1 ^ 2 /n f ) . Soli d line is 
the theoretical calcu lation from Bakes & Tielens ( 1994), dashed 
lines are that from lWeingartner & Drainel (1200 1 bb for Ry = 3.1 
(10 5 £>c = 0.0 and 6.0 for the lower and upper line), dotted lines 
are for R v = 5.5 with Case A (I0 5 b c - 0.0 and 3.0 for the lower 
and upper line; see text). 



m Horsehead 
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NGC7023 E 
tA, NGC7023 NW 
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Mon R2 

I 0.20 h 



0.15 



, £ 0.10 - 




Fig. 7. ep AH = ([Oi] 63jum + [Oi] 145 pm + [Cn] 158 pm) / 
(PAH + [Oi] 63pm + [Oi] 145 pm + [Cn] 158/fln) versus the 
charging parameter (y). The black curve depicts our calculations 
using 25 different size distributions from Weingartner & Draine 
(1200 lal) . which looks degenerated (see text). 



Our definition of e pe assumes a spherical geometry, where 
both line and continuum emissions radiate isotropically. If one 
assumes an edge-on geometry, where the continuum emission 
can escape through the molecular side while the line emissions, 
especially [0 1] 63 pm, becomes completely optically thick and 
ca n be emitted only from the front side, e pe should be divided by 
2 (lTielensll2005h . 



4.2. Comparison of e pe with theory 

Theoretically, e pe can be expressed as a function of the charg- 
ing parameter, defined as y — GqT 1 ^ 2 /^, where T is the 
gas temperature and n e is the electron density, y is propor- 
tional to the ratio of the ionization and recombination rate 
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Fig. 8. Fraction of PAH + (/(+)) versus the charging parameter 

(y). 



dBakes&Tierenslll994h . We calculated y from Go, «h> an d T 
listed in Table [4] They were carefully estimated in each subre- 
gion ( Appendix IB! to be independent of the existing PDR mod- 
els. We estimated n e by assuming that most electrons are pro- 
vided by carbon ionization and all carbon atoms are singly ion- 
ized, i.e., n e - x(C)rin, where x(C) = 1.6 x 10~ 4 (Sofi a et alJ 
2004) is the elemental carbon abundance. In Fig. [6] we show 
the relation between e pe and y. There is a trend that e pe de- 
creases when y increases. iBakes & Tielensl (1 19941) theoretically 
estimated e pe based on the physics of the p hotoelectric effect 
on dust grains. Weingart ner & Drainel(l2001bl) modeled e pe with 
some improved physical parameters and various size distribu- 
tions of dust grains in Weingartn er & Drainel (1200 lab . Dashed 
and solid lines in Fig.[6]show some of their estimates, represent- 
ing the strongest variation of e pe depending on the dust size dis- 
tribution. Dense regions characterized by R v = 5.5 (dotted lines) 
have a lower heating efficiency. The parameter be is the total 
C abundance in the log-normal populations of the size distribu- 
tion, and the model with a larger be shows a stronger enhance- 
ment at the size of < 10A dWeingartner & Drainell2001al) . Case 
A is the case without the co nstraint on the total grain volume 
dWeingartner & Drainel2001aT) . Although the observed trend that 
e pe decreases when y increases is reproduced by models, the 
observation indicates a weaker dependence of e pe on y, and all 
models overestimate e pe at low y. In the following, we focus on 
the photoelectric heating only on PAHs. 

When we consider e pe on PAHs, ([O i] 63 /mi + [O i] 145 fim 
+ [Cn] 158//m) / TIR is not an appropriate definition, because 
the TIR expresses the total energy emitted by all dust grains. 
The total energy that PAHs absorb is converted into either MIR 
AIB emission or cooling line emission through the photoelec- 
tric effect on PAHs. A fraction of the energy (generally taken to 
be 0.5) remains behi nd as electron ic excitation energy after the 
photoelectric effect dTielensll2005l) . which also results in MIR 
AIB emission. Therefore €pah is defined as ([Oi] 63 /mi + [Oi] 
145 /mi + [Cn] 158/mi) / (PAH band emission + [Oi] 63 /mi + 
[O i] 145 fim + [C n] 158 /mi). The trend of £pah against y (Fig.[7J 
is very similar to that reported in Fig. |6] 



Following iTielensi d2005l) . we computed a simple theoreti- 
cal estimate of 6pah- For a given PAH containing the number of 
carbon atoms of A^, 



1 Ihv - IP 

epAHjv f = ~f(Z = 0) — 

2 \ hv 



(1) 



where /(Z = 0) is the neutral fraction and IP is the ionization 
potential. With a typical photon energy of 10 eV and an ioniza- 
tion potential of 7 eV, this equation becomes 



ePAH,* c = 0.15/(Z = 0). 

The neutral fraction is given by 

/(z = o) = (i+ ro r 1 , 



(2) 



(3) 



where yo is the ratio of the ionization rate over the recombination 
rate, which for small PAHs is given by 



7o 



3.5 x 10^W, I/2 



(4) 



We substitute Eq.(0]i into Eq.fO and th en into Eq.flS, and in- 
tegrate over the size distribution from Weingartner & Draine 
d2001al) up to N c = 100. In contrast to e pe , in which a different 
size distribution gives a large difference (Hig. [6j, 6pah is insensi- 
tive to the adopted size distrib ution. In F ig. [7] the res ults wi th 25 
different size distributions in lWeingartner & Drainel (1200 lab are 
plotted, but they look degenerated. Although the match between 
the theoretical estimates and the observed values is not excel- 
lent, it is much better than in Fig. [6] This result suggests that the 
photoelectric heating is dominated by PAHs. 

In Fig. [8] the relation between the fraction of PAH + (/(+)) 
and y is shown. PAHs are almost neutral at y < 10 3 and almost 
fully ionized at y > 10 6 , and there is a transition in-between. 
Since y is defined by the environment properties and it is pro- 
portional to the ratio of the grain ionization and recombination 
rate, Fig. [8] confirms that the fraction of PAH + is a good indi- 
cator for the positive charging of grains in general. This trend 
is also consistent with the correlation between the intensity ratio 
of PAH (6.2 //m)/PAH(l 1 .3/mi) and y presented in Galliano et al. 
(12008b . 

5. Summary 

We analyzed HerschelfPACS and Spitzer/IRS spectroscopic ob- 
servations in six PDRs and showed that the photoelectric heating 
efficiency (e pe ) is lower in regions with a large fraction of posi- 
tively ionized PAHs (/(+)). Based on examining the photoelec- 
tric heating efficiency on PAHs, we found a dominant contribu- 
tion of PAHs to the photoelectric heating. 

Our PDR sample covers a wide range of physical conditions 
(100 < Go < 10 5 ) and provides a good test case for investigat- 
ing the relation between e pe and the charge state of PAHs. We 
estimated e pe as ([Oi] 63 fim + [Oi] 145 fim + [Cn] 158 fim) / 
TIR. e pe varies in our PDRs between 0.1 and 0.9%. /(+) was 
obtained from a fit of the MIR spectra with a set of template 
spectra representing PAH-related species and varied from 0% 
(+11%) to 87% (+10%). All positions with a high /(+) show a 
low e pe , and all positions with a high e pe show a low /(+). This 
trend supports a scenario in which a positive grain charge re- 
sults in a decreased heating efficiency. The theoretical estimate 
of e pe shows a stronger dependence on the charging parameter 
(y) than the observed e pe reported in this study, and overesti- 
mates e pe at low y. The photoelectric heating efficiency on PAHs, 
e P AH= ([Oi] 63 fim + [Oi] 145 fim + [C n] 158 fim) / (PAH band 
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emission + [Oi] 63 /mi + [Oi] 145 jt/m + [Cn] 158jum), shows 
a much better match between the observations and the theoreti- 
cal estimates, indicating a dominant contribution of PAHs on the 
photoelectric heating. PDR models that fully account for the rel- 
ative contribution of different PAH and eVSGs populations are 
needed. Velocity-resolved observations o f [O i] 63 nm in the fu- 
ture with for instance SOFIA/upGREAT (iHevmincket al.ll2012l) 
will enable us to investigate the f oreground absorption and the 
optical depth effect in detail (c.f. lBoreiko & Betdll996b . 
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beam size = 9.0" beam size = 1 1 .5" 




(2.5,0), and calculate the mean flux ratio of the adjacent spaxels 
as 



2.5 0.5 1.0 1.5 2.0 2.5 
Position of the point source 

Fig.A.l. Black curves show the simulated r, defined as 
Eq. (lA.U , as a function of the position of the point source for the 
beam size of (a) 9" and (b) 11.5". The relative contribution of 
the extended source compared to the height of the point source at 
63 //m is 0.0, 0.02, 0.1, 0.3, and 0.7 from upper to lower curves. 
Shaded areas represent the observed values with errors of ±3<x 
(blue lines in Fig. lA.2l >. (c) The flux ratio of 9" beam and 1 1.5" 
beam at spaxel 2 as a function of the position of the point source 
with the same model. Shaded areas show ranges when r in (a) or 
(b) matches the observed value. 




0.5 1.0 1.5 2.0 2.5 

Averaged flux ratio of adjacent spaxels (r) 



3.0 



Fig. A.2. Histograms of the observed r for [O i] 63 //m and [C n] 
158 /mi. Red lines show the fit with a Gaussian profile, and blue 
lines show the center and +3cr of the fitted Gaussian. 



Appendix A: Relative flux uncertainty between 
63//m and 158//m caused by the difference of 
the beam size 

We estimate the uncertainty of the relative flux in one spaxel at 
different wavelengths that arises from the different PSF and the 
unknown spatial structure of the sources by considering the com- 
bination of two extremes, pure point sources and flat extended 
emission, and by comparing the contrast between neighboring 
spaxels for these combinations with the contrast actually mea- 
sured in the observations. 

We consider three spaxels 1-3 along a line whose centers 
are located at (0.5,0), (1.5,0), and (2.5,0) in units of a spaxel 
size (9.4"). We put a point source represented by a Gaussian 
profile with the FWHM of 9" and 1 1 .5" for 63 /im and 158 /mi, 
respectively, centered at various positions between (0.5,0) and 



1 / F(spaxel 2) F(spaxel 2) 
r — — 1 

2 \F( spaxel 1) F( spaxel 3) 



(A.l) 



where F is the flux falling in one spaxel. r is lowest when the 
center of the point source is located at (1.5,0), which is the cen- 
ter of spaxel 2, and the minimum r is 7.1 and 4.0 for 63 /mi 
and 158/mi, respectively (upper lines in Figs. lA.lfa .b). On the 
other hand, we derive the observed r of the [O i] 63 //m and [C n] 
158 //m line intensities from the original 5x5 PACS spaxels for 
all regions. The histogram of the distribution of r is well repre- 
sented by a Gaussian with a small number of outliers at high r 
values (Fig. IA.2k and the fitted Gaussian has a center of 1.1 for 
both lines and cr of 0.2 and 0.1 for 63 //m and 158 /mi. Even the 
upper 3cr (1.7 and 1.4) is smaller compared to the above model 
with a point source, indicating that the contribution of extended 
emission is significant in our regions. Therefore, we modify the 
model by adding a constant extended component to all spaxels, 
and estimate r, which is shown as black lines in Figs. IA.lh and 
b with varying constant values. Then we estimate F(spaxel 2, 
63 //m)/F( spaxel 2, 158 /mi), which expresses the uncertainty of 
the relative flux at 63 //m and 158//m, for each model (black 
lines in Fig. IA.lb ). The observed r with 3cr errors are shown as 
shaded areas in Figs. IA. lb and b, and the corresponding ranges 
are shaded in Fig. IA.lt ;. which is numerically calculated by sur- 
veying the constant emission value and determining the position 
range of the point source, where r is in the observed range for 
each constant emission value. It reaches to within 15% of unity. 
Therefore, we take 15% as the uncertainty of the relative flux 
between 63 /mi and 158 /mi caused by the beam size difference. 



Appendix B: Estimate of the charging parameter 

Here we describe the estimate of Go, the gas density «h, and the 
gas temperature T to be used for deriving the charging parameter 
y in individual subregions (Table @). The rotational temperature 
derived from the low-/ H2 pure rotational emissions up to S(4) 
is adopted as T in all regions. Go and «e affect y linearly and 
we estimate them independent of an existing PDR model. These 
properties are representative of the warm PDR surface where 
most of the H2 rotational emission and PAH emission comes 
from. 



B.1. Horsehead 



Habart et al.l (120051) estimated G (star) ~ 100 at the PDR inter- 
face based on the effective temperature of the exciting star and 
assuming geometrical dilution at the projected distance. Under 
this assumption, our three subregions are close enough to each 
other to make no difference for the Go(star) between them. Go 
derived from TIR, G (TIR), is listed in Table H We adopt both 
Go(star) and Go(TIR) for the possible range in individual sub- 
regions. For the gas temperature T, we use the range of the es- 
timate from the intensity ratio of H 2 S(2)/S(0), S(3)/S(l), and 
S( 4VS(2) in|Habart et all l20Tlh . 

lAbergel et al.l (120031) showed that the lower limit of the den- 
sity behind the fil ament is ~ 2 x 10 4 cm -3 from the infrared 
brightness profile. Hab art et al.l (|2005) modeled the spatial dis- 
tribution of H2, PAH, CO and 1 mm continuum emissions us- 
ing a PDR model and suggested the density gradient from «h= 
10 4 ctrT 3 in the H2 emitting region to «h= 2 x 10 5 cirT 3 in the 
inner cold molecular layers. Since the region we are interested 
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in is the H 2 -emitting gas, we still include the «h value in the 
H2-emitting region even for subregion 3, i.e. we adopt «h=(0.1- 
2)xl0 5 cm 4 for subregions 2 and 3. «h~ 10 4 cm" 3 at the inter- 
face is also confirmed by the pressure equilibrium with the ion- 
ized gas. In the ioni zed gas, n e is estimated to be 100-350 cm -3 
with T e = 7500 K (Compiegne et al. 2007). Assuming the pres- 
sure equilibrium, 2n e T e -n\(I ', gives 71h=(0.3-2)x10 4 cm" 3 . 
Therefore, we adopt «h=(0.3-2)x10 4 cm" 3 for subregion 1. 

B.2. Ced201 

iKemperet all d 19991) estimated Go = 200 by comparing the ob- 
servations with a PDR model in units of the average interstel- 
lar radiation between 2 eV and 13.6 eV. On the other hand, 
lYoung Owl et al.l d2002) derived Go = 300 by the infrared con- 
tinuum emission, with a geometry correction factor of 2. We 
adopt Go = 200 as a lower limit and Go (FIR) as an upper 
limit for each subregion. The fraction of eVSG s determined by 
PAHTAT is correlated to G (iPilleri et al.l2012bl) and the adopted 
Go for each subregion is also consistent with it. 

For n H , lYoung Owl et al.1 (|2002|) derived n H = 4 x 10 2 cm" 3 
from t he lower limit of [C n] 158/im/[Oi] 63 /im using a PDR 
model. Kempe retal] (1 1 9991) estimated n H of (5 + l)xl0 3 cm" 3 by 
simple excitation models of CO and I3 CO emissions. They also 
modeled the emission of CO, C, C\ CS, and HCO + using a PDR 
model, which suggests «h of 1.2 x 10 4 cm" 3 . We use all these 
ranges for the two subregions. T is esti mated as ~ 330 K from 
the ratio of H 2 S(l) and S(3) emissions (iKemper et al.ll 19991) . 

B.3. NGC 7023 E 

IPilleri et all (1201 2bl) modeled the spatial distribution of the MIR 
AIB emission, assuming a spherical shell geometry and the en- 
ergy balance taking into account, and derived the spatial varia- 
tion of n H and G along a cut in NGC 7023 E. At the PDR front, 
Go is calculated to be 250 based on the spectral type of the star, 
geometrical dilut ion and assuming a n attenuation of A v = 1.5 
around the star (IPilleri et al.ll2012bl) . The density at the PDR 
front is derived as 1 .4x 10 3 cm" 3 . Around the MIR AIB emission 
peak, which corresponds to the regions in this study, Go is esti- 
mated to be 120-170 and n H = (0.5-1) x 10 4 cm" 3 . On the other 
hand, G (TIR) is 970 and 900 for two subregions. We adopt 120 
as a lower limit of Go, and the corresponding Go(TIR) as an up- 
per limit for each subregion, and «h= (0.1-l)xl0 4 cm" 3 for both 
subregions. For the gas temperature T, we use 258 -370 K from 
the in tensity ratio of H 2 S(2)/S(0) and S(3)/S(l) dHabart et all 
1201 lh . 

B.4. NGC 7023 NW 

The sa me modeling as for NGC 7023 E was made bv lPilleri et al.l 
(2012b). Our three subregions are not exactly on their cut be- 
cause we chose the subregions to maximize the variation of 
the P AH + fraction. Subr egion 2 is located at the PDR inter- 
face. IPilleri et al.l (l2012bl) estimate G (star)= 2600 at the PDR 
front, whereas Go(TIR) — 7700 at this subregion. We adopt 
Go = 2600-7700 for subregion 2. Subregion 1 is located closer 
to the star than the PDR interface. The scaling of Go(star) by the 
projected distance gives Go = 10 4 , higher than Go(TIR), pos- 
sibly because either the real distance to the star is larger than 
the projected distance, and/or the assumption that all UV radi- 
ation is absorbed and converted into the IR emission underesti- 
mates Go in such a cavity because many UV photons pass the 



region. Nevertheless, we conservatively cover both values, i.e., 
Go = (0.5-1) x 10 4 , for subregion 1. The adopted Go in sub- 
regions 1 and 2 is also consistent wit h the correlation bet ween 
the fraction of eVSGs and G from IPilleri et al.l (l2012bl) . For 
subregion 3, we us e this correlation to derive Go (see Eq.(5) in 
IPilleri et al.ll2012bl) . which results in G = 500-3100. Together 
with G (TIR)= 5400, we a dopt G = 500-5400 for subregion 3. 

In lPiiieri et al.l (l2012bl) . the modeled density profile sharply 
increases from 1.1 X 10 3 cm" 3 to 2 x 10 4 cm" 3 within ~ 8". 
iFuente et al.l {1996) showed high-density filaments of a few 
10 5 cm" 3 based on HCO + observations. Therefore, we adopt 
10 3 -2 x 10 s cm" 3 for subreg i on 2 and (0.2-2)xl0 5 cm" 3 for 
subregion 3. Berne & Tielen sf (120 121) examined nu in the cavity 
in detail from several diagnostics and suggested 150+100 cm -3 . 
We take this as «h for our subregion 1 . For the gas temperature 
T, we use 43 0^450 K from the intensity ratio of H 2 S(3)/S(l) 
and S(4)/S(2) (IFuente et alJl2000h . 

B.5. Carina N 

We estimate Go(star) by compu ting the contr ibution from all OB 
stars of Trumpler 14 listed in ISmithl (120061) . which gives (7- 
8)xl0 3 for the three subregions. On the other hand, Go(TIR) 
is (1.5— 1.6)xl0 . The point-sourc e catalog of the W ide Field 
Infrared Survey Explorer (WISE; IWright et al.ll2010b contains 
several sources with a flux stronger at 4.6 yum than at 3.4/mi 
close to our region, which indicates a contribution from embed- 
ded protostars. We consider both Go(star) and Go(TIR) as a pos- 
sible range in individual subregions. 

Kram er et a D d2008) modeled CO and [C i] emissions in 
the Carina N region and suggested «h of 2 x 10 5 cm" 3 . This 
is consistent with a pressure equilibrium with the H n region. 
Using the emission lines ratio from ions with similar ioniza- 
tion potential observed with IRS, we can estimate n e assum- 
ing the elemental solar abundance and isothermal thin emission. 
[Nen] 12.8 fxm/[S m] 18.7 /im, [Pra] 17.9//m/[Nen] 12.8/mi, 
and [Arm] 8.99//m/[Sm] 18.7 /xm give n e of 10 4 , 2 x 10 4 , and 
9 x 10 3 -10 4 cm" 3 , respectively, when T e = 10 4 K. Since the gas 
temperature T in PDRs is derived as T — 290-529 K from H 2 
emission lines, «h can be estimated as 3 x 10 s — 10 6 cm" 3 . We use 
2 x 10 s — 10 6 cm" 3 for all subregions. 

B.6. Mon R2 

We calculate Go(star) from the contribut ion of IRS 1-4 (see 
Fig. Q]) with the luminosity listed in iHenning et alJ dl992l) . as 
(5.2-5.9)xl0 4 , 10 5 , and (l.l-1.2)xl0 5 in the three subregions. 
While Go = 5 x 10 5 i s used to charac terize the radiation field 
at the ionization front dRizzo et al 1 120031) . our subregions do not 
cover the nearest regions of infrared sources to avoid unresolved 
complex spatial distributions. Go(star) matches Go(TIR) well. 
We consider both Go(star) and Go(TIR) as a possible range in 
individual subregions. 

iRizzo et alJ (120051) and IPilleri et al l (120131) performed an 
LVG analysis with C3H 2 an d C?H at several posi tions in Mon R2 
and derived «h > 10 6 cm" 3 . iGinard et all (120 1 2|) showed similar 
results from different molecules; the molecular hydrogen den- 
sity of a few 10 5 cm" 3 to 10 6 cm" 3 . These molecules traces 
the cold molecular gas wit h a temperature of T ~ 45 K 
(iGian nakopoul ouet al.l 1 1 997h . On the other hand. iBerne et al.l 
(2009]) derived w H of (0.4-4)xl0 5 from PDR model calculations 
using H 2 rotational emission lines, and the temperature derived 
from the H 2 rotational emission lines is T — 314-574 K. These 
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«h and T are consistent with a pressure equilibrium with the 
cold molecular gas. In the Hn region, the emission measure of 
1.9 x 10 7 pccm~ 6 and the geometrical mean diam eter of 23.7" 
dTakahashi et al.ll2000t IWood & Churchwellll 19891) give an esti- 
mate of n e = 1.4 x 10 4 cm" 3 . With T e = 7600 K dDownes et alj 
fl97l and T — 314-574 K, the pressure equilibrium gives «h 
=(4-7) xlO 5 cm -3 , which is also compatible with the estimate 
from H2 rotational emissions, although the pressure equilibrium 
is an inadequate assumption in an UCH 11 like Mon R2. We adopt 
the estimate range from H2 emissions, i.e., (0.4^4)xl0 5 for all 
subregions. 
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